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SUMMARY 


A technique was formulated for constructing spectral reflectance estimates 
from fflultispectral data obtained with the Viking lander cameras . The output 
of each channel was expressed as a linear function of 'the unknown spectral 
reflectance producing a set of linear equations which were used to determine 
the coefficients in a representation of the spectral reflectance estimate as 
a natural cubic spline. The technique was used to produce spectral reflectance 
estimates for a variety of actual and hypothetical spectral reflectances. 


INTRODUCTION 


The Viking lander cameras (ref. l) will return multispectral images of 
the Martian surface with four orders of magnitude higher resolution than has 
been previously obtained. It is desired to extract spectral reflectance curves 
from this data. However, the data are limited to 6 spectral channels and most 
of these channels exhibit out-of-band response. 

It is inappropriate to generate a data point for each channel by 
associating a reflectance value with a distinct wavelength; this is particu- 
larly true for those channels with appreciable out-of-band response. It is 
•unlikely that data points so constructed will lie on the true spectral 
reflectance curve, -and that any method of fitting a curve to these points 
will adequately approximate the true reflectance. 

Instead the output of each channel can be expressed as a linear integral 
function of the unknown spectral reflectance and the known solar irradiance, 
atmospheric transmittance, camera optical throughput, and channel responsivity. 
This produces 6 equations - one per channel - which can be used to determine 
the coefficients in a representation of the spectral reflectance as a natural 
cubic spline. In this paper the appropriateness of this technique is demon- 
strated by using it to produce accurate approximations to the true spectral 
reflectance of 8 materials felt likely to be present on the Martian surface 
and l6 hypothetical spectral reflectances chosen for illustrative purposes. 
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FOEMULATIOH 


Let p(X) denote the (nnknown) spectral reflectance at -wavelength X of 
the material that is imaged by the Viking lander camera. Knowledge of p(A) is 
limited to 6 spectral samples. Except for a channel-dependent, multiplicative 
constant, which can he determined by a calibration using as reference a test 
chart on board the Viking lander (see ref. 2), these 6 spectral samples are 
given by b^ where 


Li = ) T^(x)p(x)dX 


i = 1,2, . . .,6 (l) 


The system transfer f-unctions T^(x) are given by 


S(X) t(x) rAx) R.(X) 
tAx) = ^ ^ i 


i = 1,2, .. .,6 


where S(x) is the solar irradiance, Ta(x) the atmospheric transmittance, 
T(,(X) the camera optical throughput, Rj_(x) the channel responsivity, and t^ 
is a constant chosen so that 



i = 1,2, .. .,6 


Plots of typical system transfer functions are shown in figure 1. Note 
specifically the appreciable out-of-band response of the Blue (i=l), IE2 
(i=5)j and IR3(i=6) channels. Note also that with the possible exception of 
the Green (i=2) channel, none of the system transfer functions are adequately 
approximated by an impulse f tine t ion. 


THE REFLECTANCE ESTIMATE AS 
A NATURAL CUBIC SPLINE 


Equations (l) describe the relationship between the 6 multispectral 
samples bp, b2, ...,b5 and the unkno-wn reflectance p(x). These six equations 
can be used to produce a natural cubic spline estimate of p(x), denoted 
< p (x ) > j where 


7 

< p(x) > = X. C(X-X.) (2) 

j=0 ^ 
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The 8 Imots \q, are chosen to he equally spaced and located at 

the wavelengths 


X. = .33 + jA 
J 


2 ~ 0 ? 2 , 


• • } 
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where the spacing is A= .12 |j,m. Recall that each cubic spline basis function 
C(l-\j) is a bell-shaped curve centered at the knot and defined by C(x) 
where ^ 



The coefficients xq, Xj_^...^xy are to be determined. 

It is desirable to impose the natural boundary conditions < p(x) >" = 0 
at the knots Ip and Xg. These two conditions give rise to the equations 


Xq - 2x^ + x^ == 0 (3a) 

and 

V 

^5 - + x^ = 0 (3b) 

The remaining six equations which determine the 8 coefficients are obtained by 
requiring the estimate < p (x ) > and actual reflectance p (X ) to have in- 
distinguishable camera multi spectral responses (ref. 2), i.e.. 


J“T^(x)<p(x)>dX = J”T^(x)p(x)dX 


l_y 2^ . • . ^ 8 




This produces the six equations 
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• • • ^ 


(5) 
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a. . 


X. = To. 
3 1 


i = 1,2, 


where is given by equation (l) and 




T. (X) C(X - X ) dX 


(6) 


To summarize, a natural cubic spline reflectance estimate corresponding 
to the multispeetral sample b^, b2, . . . ,bg can be produced as follows : 


(i) 

(ii) 


evaluate the 6x8 = 48 coefficients 
determine the 8 coefficients x . 


(3b), and (5) 


a. . given by equation (6) 
by solving equations (3a), 


(iii) form the estimate <p (x)> given by equation (2) 


The estimate constructed in this manner reduces to an interpolating spline 
in the idealized situation where each system transfer function can be repre- 
sented as an impulse function. To see this suppose that 


T^(X) = 5(X -X^) i = l,2, ...,6 

where the impulse system transfer functions occur at the discrete wavelengths 
Xi^ In this special case 


and 


b^ = I 6(X -X^) p(x) dX = p(X^) i = l,2,..,,6 

•'o 


a. . 



-X. ) C(x -xjdx 

-*• J 


C{Xi-X.) 


so that <p(x)>' is the (unique) nattiral cubic spline which interpolates the 
spectral reflectances p (X^)^ p (Xg), • . p (Xg) • 
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RESULTS 


Reflectance estimates were computed for 8 materials felt likely to Le 
present on the Martian surface and for l6 hypothetical spectral reflectances 
chosen for testing and illustrative purposes. These estimates are presented 
in figures 2, 3, and it. In each case the actual spectral reflectance is shown 
as a sequence of T1 discrete points (circles) in the wavelength range 
.h < X < 1.1 |j,m and the corresponding estimate is shown as a continuous curve. 
For each of the spectral reflectanc.es the corresponding multispectral sample 
bi,h2, ...^hg was calculated from equation (l) using a 71 point Simpson's Rule. 
The coefficients aj_j were calculated in the same manner from equation (6) 
assuming each system transfer function to he zero outside the effective range 
of the camera photosensor arrays^ namely < X < 1.1 pm. 

Figure 2 illustrates the reflectance estimates, for the 8 materials felt 
likely to he present on the Martian surface. For those 5 simple reflectances 
(i.e., plnacetes 5 aJid 2QA, Syrtis Major, augite, and average Mars) the 
estimates are excellent. For the 3 more complex reflectances (i.e., limonite, 
hypersthene, and olivine) the estimates are very good. The dominant features 
are reproduced^ however, due to the undersampling inherent with just 6 
channels, small period features are lost. Note particularly that the dominant 
absorption hand (at \ .93 pin) for hypersthene is quite accurately estimated. 

Figure 3 illustrates the reflectance estimates for 8 hypothetical spectral 
reflectances. The first k of these spectral reflectances (3a, 3h, 3c, and 3d) 
are very smooth and the corresponding estimates are almost exact. The next 
four (3e, 3f, 3g, and 3h) are not smooth hut the estimates remain good. Note 
specifically that the pronounced minima in 3e 8-nd 3t accurately reproduced. 
Note also in 3g the characteristic oscillation exhibited by the natural cubic 
spline in the neighborhood of a large slope. In 5h the loss of small period 
features is again evident. 

Figure 4 illustrates the reflectance estimates for 8 additional hypo- 
thetical spectral reflectances . All 8 of these are of the form 


p (x) 


.25 + 



(7) 


where the parameters a, p have the values : 



figure 4 
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As the period becomes shorter (i.e., as p decreases), the quality of the 
reflectance estimates deteriorates. This is particularly evident in the 
sequence h-B., kc, 4e, 4g and less evident in the sequence 4b, 4d, 4f, 4h. It 
is also true that the quality of the estimate is affected by the location of 
dominant spectral reflectance features relative to the location of the system 
transfer functions. This is illustrated by figures 4e and 4f where the 
spectral reflectances differ only by a shift of .07 (j,m while the corre- 
sponding estimates differ dramatically. Figure 4g is a clear demonstration 
of aliasing whereby a short period harmonic spectral reflectance curve has a 
reflectance estimate which is nearly harmonic but with a larger (false) period. 


CONCLUDING REMARKS 


A technique was formulated for constructing natural cubic spline spectral 
reflectance estimates from multispectral data obtained with the Viking . lander 
camera. Using this technique it was demonstrated that smooth, simple spectral 
reflectance curves can be estimated almost exactly. For more complex spectral 
reflectance curves, large period features can be faithfully reproduced; small 
period features are lost due to the undersampling inherent with the limited 
number of spectral channels . The technique completely compensates for system 
transfer functions with irregular shapes and appreciable out-of-band trans- 
mittance. Moreover the technique should be a valuable aid in selecting the 
number of spectral channels and their responsivity shapes when designing a 
multispectral imaging system. This design approach would prove to be of value 
especially if spectral reflectance properties of interest are known a priori 
and if the transfer function shapes are desired to be broad to obtain good 
signal-to-noise ratios. 
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